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( — ^ ■ Abstract 

In this work we study the electromagnetic field at Finite Temperature via the massless DKP for- 
malism. The constraint analysis is performed and the partition function for the theory is constructed 
. and computed. When it is specialized to the spin 1 sector we obtain the well-known result for the 

D ■ thermodynamic equilibrium of the electromagnetic field. 

: 

1 Introduction 

J> ■ Quantum Field Theory at Finite Temperature was motivated by the increasing interest in studying the 

properties of matter under extreme conditions as, for example, at very high temperature or density. 
The pioneering works joining together the Statistical and Quantum Field Theory were developed mainly 
by Matsubara pQ in a non relativistic context and, the relativistic case by Fradkin [2] who, via the 
functional approach, studied the different methods for calculating the Thermal Green's functions as well 
' as the structure of the Ward identities in QED4. Later other works within the Thermal Field Theory 

appeared [3] whose principal interest was to explore the possibility of restoring some broken symmetries 
that occur at zero temperature, as for example the SU(2) x U(l) symmetry of weak interaction. The 
Finite Temperature gauge theories and the problems concerning to the choice of a physical gauge and its 
O 1' dependence was analyzed by Bernard 0], in particular, the free electromagnetic field. 

On the other hand, at zero temperature, there is an alternative way to study the properties for the 
electromagnetic field which is known as the massless Duffin-Kemmer-Petiau (DKP) theory [5] that is 
not a trivial limit of the massive Dufhn-Kemmer-Petiau theory (DKP) [6] that appears as an alternative 
formalism for the description of the spin and spin I particles in a unified formulation. The DKP theory 
gives a first order linear equation and it is very similar to the Dirac one but the matrices satisfy a 
different algebraic relation. The massive and massless case of the theory were considered in the works 
[7J [HI [S] where the equivalence of the DKP theory with the theories like Klein-Gordon-Fock (KGF) and 
Maxwell was proved in Minkokwski space-time [7j and studied in curved space-time such as Riemann [5] 
and Riemann-Cartan [5]. 

At Finite Temperature the massive case is treated in the work [10] where the Bose-Einstein condensa- 
tion is investigated in the spin sector. The equivalence of many-photons Thermal Green's functions of 
the DKP and KGF theories was also proved for the scalar sector [11] calculating the polarization operator 
at 1-loop order and, in |12j is shown the equivalence of many-gluons Green's functions in the DKP and 
KGF Statistical Quantum Field Theories 

As above mentioned, all accomplished studies on the massless DKP theory were made at zero tem- 
perature. The aim of this work is to study the thermodynamics of the electromagnetic field by using the 
massless DKP theory. The paper is organized as following: In section 2, we give a brief review of the 
massless DKP formalism considering a theory with only one real DKP field. In section 3, the constraint 
analysis is performed for the DKP theory and it has been shown that the model has two first class 
constraints. In section 4, we calculate the partition function and, finally, we give our conclusions and 
commentaries. 
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2 The Massless DKP Theory 

The massless DKP theory is described in Minkowski space-time by the following Lagrangian density [5] 

C M = i$ j^d^ - - ^ , (1) 

where ?/> = ip^rjo with 770 = (2/3q — l). The /3 M and 7 are singular square matrices satisfying the following 
algebra 

/F/3 A /T + v /3 x f3>* = t/^/T + , (2) 

Fl + lP l = P t , 7 2 =7 ■ (3) 

Due to the singular character of the /3 M matrices the transition of the massive case to the massless theory 
is non trivial and demands a different treatment. On the other hand the representations of the DKP 
algebra are reducible and contain the sectors of spin and spin 1 in their structure. 

From the Lagrangian (fTJ) we obtain the equation of motion for the massless DKP field 

1/3^8^-^ = . (4) 

It can be shown that the Lagrangian |T]) and the massless DKP equation remain invariants under the 
following gauge transformation 

tl> -» V' = V> + (1 - 7) $ 1 (5) 

iff the field $ satisfies the condition 

i^fy (1 - 7) $ = . (6) 

When the fields under consideration are no charged we have a real DKP fielcfl ip, in such situation 
the Lagrangian ([TJ takes the following form 

C = \i? drf - % -d^ T (7/^7) i> - (ry° 7 ) i, (7) 

3 Constraint Analysis 

We proceed the study of the constraint analysis to the real massless DKP theory from the Lagrangian 
([7]) which is written as 

c = l - r (v°i^) ab d^ b - \d^ a (r,°(3^) ab ^ - l -r ( v ° 7 ) ab ^ b . % 

as usual we define the canonical momentum ir a as 
from which a set of primary constraints appear 6 

a = 7r a +i(/?°7) oi y , (io) 

because the two different representations for the (3^ matrices we have for the spin sector that a = 
ip, 0, 1, 2, 3 and for spin 1 sector that a = 0, 1, 2, 9. 

The canonical Hamiltonian density TLc that follows from the Lagrangian ((5J) is given by 

n c = \d k r {v°P k i) ab ^ b - \r ( W fe ) Q& d k ^ b + \r (v°i) ab v b (u) 

and considering the set of constraints (|10[) we have the primary Hamiltonian density Tip as 

H P =Hc + \ a e a , (12) 

where A a are the Lagrange multiplier. 

The Poisson bracket (PB) for the primary constraints results in 

{MxMb(y)} = *(/3°W(x-y) • (13) 



1 When the DKP field is real, the (3^ matrices must been satisfied 0q = — /3rj , 0£ = flu ■ The representation of the (3^ 
matrices for the spin 1 sector is given in the appendix A, where we also included the 7 matrix. 
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To investigate the possibility of obtaining more constraints in the theory we apply the preservation in 
time of 6 a , i.e. 

o (x) = {0 o (x),tf P } , (14) 



where H P = ]d.H P is the primary Hamiltonian. Thus, the stability condition provide 

9 a = i(r, o /3 k ) ab d k i; b -( V o 1 ) ab Tp b + o ab \ b ttO , (15) 

being (3^ singular matrices we conclude that not all A coefficients can be obtained from the relation (fTS"]) 
and more constraints appear. These new constraints are selected by means of the projector M = (l — , 
thus, we obtain 

6^ = (M(3 k ) ac zd k r - (M 7 ) ac V c (16) 

that are a set of secondary constrains. When the preservation in time of these secondary constraints is 
imposed no more constraints appear in the theory 

Now we calculate the PB for all primary and secondary constraints 

{MxMf(y)} = [(i(3 k d% + 7 )M] afc( 5(x- y ) , 

{0f(x)A(y)} = [M(z/3^- 7 )] a ^(x-y) , (17) 



{O^ (xM£ 2) (y)} 







We can write all set of constraints as C a — \ a , 6 a > such that its matrix is 



{Ca(x),C & (y)} = 



i0° (ip k d%+j)M 
M(i/3 k d£ - 7) 



*(x-y) (18) 



and the determinant of {£ a , C,b\ is zero. But the constraint structure depends of the chosen representation 
for the DKP algebra as we will see go on. 



3.1 Spin 1 sector! 



For this sector the representation of the /3 M matrices is 10 x 10 and the spin 1 DKP field is a column 
matrix with ten real components 

if>= (^\^^^^V^^^ 9 ) T ■ (19) 

Consequently, in this sector, from (|10p we observe that there are ten primary constraints and write 
explicitly 

0a = ir o , (20) 

01=7n+^ 7 , 6 2 =TT 2 +ij 8 , 03=^3+^ 9 , (21) 

e n = ir n , n = 4,5,6,7,8,9 . (22) 
And from (|16|i we obtain four secondary constraints given by 

0<j 2) = -d^ 7 - d 2 ip & - d 3 i/j 9 , (23) 



q(2) 



- d 3 ip 2 - 9 2 V 3 - V' 4 , = diip 3 - dsip 1 - V' 5 , 0f ] = d 2 ^ - dnp 2 - V- 6 , (24) 



To classify these constraints as first and second class we perform the calculation of the PB between 
all these primary and secondary constraints such as it is shown by the matrix {( a (x) , (y)} in (fl"g|) . 
In our case, the rank of the matrix {Ca,Cb} is 12 which is the number of second class constraints and, 



2 For the spin 0, the representation for the /3 M matrices is 5 X 5. The rank of the constraint matrix {£ a , is 8, the null 
space has two trivial constraints which are irrelevant, and the set of constraints is second class. 
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the dimension of the null space is 2 that gives the number of the first class constraints. The null 
space is formed by the constraint 6q = 7i"o and by the linear combination of second class constraints 

(2) 

di9i + c>2#2 + 936*3 + #o that defines another first class constraint G — dk^k- Thus, we obtain two first 
class constraints 



7TQ 



G = <9fc7Tfc 



1,2,3 



(25) 



and twelve second class constraints given by the equations (|2T1) . ([22)1 and ([24 
The projectors of the spin 1 sector are defined as 

R» = (/3 1 ) 2 (/3 2 ) 2 (/3 3 ) 2 

R n» = Ril p ) M ,i/ = 0,1,2,3. 



(26) 



such that the field ip^ = R^ip is a Lorentz vector and ip^ v = R^ip is an antisymmetric second-rank 
Lorentz tensor; when multiplied with the 7 matrix we also get R^jtp = and R^-ftp = R^ip. 

Using the projectors and R^ v and, from the relation ([5]) we conclude that only the vector compo- 
nents of the DKP field are transformed as it is shown to follow 



= , 



4,5,6,7,8,9 



(27) 



and from ([6]) we get = $ M = ±<9 M A, being A an arbitrary scalar function, thus, we can conclude 
that the theory under consideration is a local U (1) gauge field theory. Then, we impose the following 
gauge fixing conditions 

fii = 0kV* . ^2 = ^° , (28) 
such that the set xa> = {&q, G, Oi , f^} is second class, thus, the PB matrix of the set is given by 



Da'B' = {xa> (x) ,Xa< (y)} = 












-1 








-A* 








A x 








1 












<5(x-y) 



(29) 



where A = (dk) = (di) + (82) + {83) , computing the functional determinant we get 



det D A > b> = [det A] 



(30) 



4 The Partition Function 

Now we study the thermodynamic equilibrium of the electromagnetic field using the DKP formalism. 
Such as we see the constraint analysis of the theory gives for the spin 1 sector two first class constraints 
characterizing a local U(l) gauge theory. Then, we write the partition function for the massless real DKP 
field using the Hamiltonian formalism 

Z= J V4,VttS(Q a ) S (xa') (det C AB ) 1/2 (det D A , B ,) 1/2 exp [J d 4 x (iir a d T ^ a - H c )} , (31) 

periodic 

where He is given by ([TTj) and 

T>i]) = Vil) a , Vn = VTr a , a = 0,1, 2,..., 9 . (32) 
The fields ip are restricted by the periodicity condition 

^(0,x)=V(/3,x) , (33) 

where Q A — |#a, ^j, 2 ' j is the set of second class constraints given by the equations ([2"Tj) . ([22]) and |24[) ; 

the set XA' = {#0, G, fix, ^2} is given by the set of first class constraints ([2"S"[) and its respective gauge 
fixing conditions (|28[). The matrix C A b = {9,4,0s} can be obtained from (p~8|) and its determinant is 
deters = 1 . The matrix Da'B' = {xa 1 ,Xb'} and its determinant det D A ib> are given by the equations 
(f29[) and pO]) . respectively. 
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But, it is interesting to perform the calculation of the partition function in a manifest covariant way. 
Thus, it is possible to show that the equation (|3ip becomes 



Z = N((3) j VijjS (F [ifj A ] ) det 

periodic 



Qpg 



dA 



exp 



d A x 



(if3 A d A - 7) V> 



(34) 



where F \jp A1 \ is an arbitrary gauge fixing condition. Here we consider 



F° [^ A ]=F[^]--^DA 



(35) 



with the gauge transformation ip A — > i\) A — d A A , and / is an arbitrary scalar function. 

It is worthwhile to note that (|34| is exactly the Faddeev-Popov technique [14] used to quantize a local 
gauge theory. Consequently, the equation ([34]) can be expressed as being 



Z = N{I3) J X^det 

periodic 



□ 



exp < d x 



(36) 



where the index A — r, 1,2,3 and (3 A d A — if3°d T + (i k d k . We remark that at zero temperature the 
det (□) is a constant that can be ignored, however, at Finite Temperature it turns out a very important 
temperature dependent term. Using the projectors R A we rewrite the gauge fixing term in matrix form 
such that the partition function reads as 



Z = N{(3) det 



Z' 



where 



Z' = J 2?^ exp | j d A x i^V (iP A d A (R A f R B d A d B ^j ^ 



(37) 
(38) 



periodic 

and R A d A = iR°d T + R k d k . 

To perform the calculation of the functional integral we use the Fourier series of the DKP field 



(39) 



with Lo n — 2im/l3 and the periodicity conditions ([3"3"]) for the DKP field are imposed. 
Then, substituting ([39]) in ([38]) we get 



Z' 



J|det 



icj n p° + Pk p k + 7 + ~?7 (iu n R° + PkR k f {iu n R° + PkR k ) 



-1/2 



n 



c 



-1/2 



n 

n , p 



(40) 



where tu p = |p|. The determinant in (|37jl is 

det 

Finally, the partition function reads 



] r K + ^p) 1 

n rt V ~ 



Z = JV(/3)I]K+^)- 1 , 



n,p 



as we are interesting in InZ, then, from the last expression we obtain 



(41) 



(42) 



(43) 
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The value for the N {(3) is selected in a manner that it cancels the divergent term, i.e we choose 



next, we perform the sum (see, for example, [4l I15j) and get the following expression for the partition 
function 



which describe a massless bosonic field with two polarization states that is the characteristic behavior of 
the electromagnetic field in thermodynamical equilibrium. 

5 Conclusions 

In this work we study the massless DKP theory at Finite Temperature, and it is shown the constraint 
structure of the model leads to conclude that it is a local U(\) gauge theory in its spin 1 sector. Such 
analysis allow to construct the correct partition function using the Hamiltonian procedure. Also, we show 
that it is possible to arrive to the covariant expression which is exactly the covariant quantization of a 
gauge theory using the Faddeev-Popov approach. Consequently, the partition function of the spin 1 sector 
gives the partition function of a zero-mass Bose gas with two polarization states, i.e. the electromagnetic 
field modes in thermodynamical equilibrium. 

The perspectives to be followed are to study the Finite Temperature properties of the massless DKP 
field in curved space-time and, consequently, analyze the curvature effects in the thermodynamics of the 
electromagnetic field via the DKP formalism. Advances in this directions will be reported elsewhere. 
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A Spin 1 representation for the massless DKP algebra 
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